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ABSTRACT 

This is the second of a series of papers investigating the oscillation properties of rela- 
tivistic, non-selfgravitating tori orbiting around black holes. Extending the work done in a 
Schwarzschild background, we here consider the axisymmetric oscillations of vertically inte- 
grated tori in a Kerr spacetime. The tori are modeled with a number of different non-Keplerian 
distributions of specific angular momentum and we discuss how the oscillation properties de- 
pend on these and on the rotation of the central black hole. We first consider a local analysis 
to highlight the relations between acoustic and epicyclic oscillations in a Kerr spacetime and 
subsequently perform a global eigenmode analysis to compute the axisymmetric p modes. In 
analogy with what found in a Schwarzschild background, these modes behave as sound waves 
that are modified by rotation and are globally trapped in the torus. For constant distributions 
of specific angular momentum, the eigenfrequencies appear in a sequence 2:3:4:... which is 
essentially independent of the size of the disc and of the black hole rotation. For non-constant 
distributions of angular momentum, on the other hand, the sequence depends on the properties 
of the disc and on the spin of the black hole, becoming harmonic for sufficiently large tori. 
We also comment on how p modes could explain the high frequency quasi-periodic oscilla- 
tions observed in low-mass X-ray binaries with a black hole candidate and the properties of 
an equivalent model in Newtonian physics. 
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1 INTRODUCTION 

A recent series of papers (Font & Daigne, 2002a,b; Daigne & 
Font, 2004) has shown, through general relativistic hydrodynamic 
simulations, that the relativistic geometrically thick discs of high- 
density matter expected to form after the coalescence of a binary 
system of neutron stars or after the gravitational collapse of a ro- 
tating supermassive star may be subject to the runaway instability 
which could destroy the disc on a dynamical timescale (Abramow- 
icz, Calvani and Nobili, 1983). Accretion tori of much smaller 
rest-mass densities have also been observed in three-dimensional 
general relativistic magnetohydrodynamic simulations investigat- 
ing the properties of accretion flows onto Kerr black holes (De Vil- 
liers et al. 2003). 

With the aim of assessing whether the onset of the runaway 
instability depends sensitively on the choice of initial conditions, 
Zanotti et al. (2003) carried out general relativistic hydrodynamic 
simulations of perturbed tori orbiting around a Schwarzschild 
black hole and shown that the instability may take place indepen- 
dently of the way the mass accretion is induced. In addition, the 
simulations have indicated how the introduction of perturbations 
triggers harmonic oscillations in the tori which could produce large 



variations of their mass quadrupole. In the case in which the torus 
is made of high-density matter (as in the case of binary neutron star 
merger), these oscillations could then lead to the emission of grav- 
itational waves with amplitudes comparable with those produced 
in a gravitational stellar-core collapse, thus making these objects 
promising sources of gravitational radiation even if the instability 
does not set in. 

Together with the general relativistic hydrodynamic simula- 
tions, oscillation modes of geometrically thick relativistic discs can 
also be studied through perturbative analyses. Clearly, this second 
approach is computationally less intensive and allows therefore for 
a more detailed investigation of the parameter space. Using a linear 
perturbative approach, Rezzolla et al. (2003b; Paper I hereafter), 
have recently investigated the axisymmetric oscillation properties 
of relativistic tori in a Schwarzschild spacetime. To simplify the 
treatment and make it as analytical as possible, the tori were built 
with vertically integrated and vertically averaged quantities, thus 
transforming the eigenvalue problem into a set of coupled ordinary 
differential equations. The tori were then modeled with a number 
of different non-Keplerian distributions of specific angular momen- 
tum. Overall, the perturbative analysis of Paper I confirmed the re- 
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suits obtained by Zanotti et al. (2003) through nonlinear simula- 
tions, revealing that the oscillations found in the simulations corre- 
spond to p modes and that the lowest-order eigenfrequencies must 
be in a sequence of small integers 2:3:4..., rather independently of 
the size and the specific angular momentum distribution of the tori. 
Furthermore, a detailed study of the eigenfunctions in Paper I has 
also provided a possible explanation for the numerical evidence that 
the runaway instability can be efficiently suppressed for tori with 
distributions of specific angular momentum that follow a power- 
law (Font & Daigne, 2002b). 

This paper intends to extend the work of Paper I and is there- 
fore devoted to both a local and a global perturbative analysis of 
axisymmetric modes of oscillation of relativistic tori in the back- 
ground spacetime of a Kerr black hole. While the local analysis 
provides us with the dispersion relation for inertial-acoustic waves 
in relativistic non-Keplerian discs, the global approach provides us 
with the eigenfunctions and eigenfrequencies of the p-mode oscil- 
lations of the system. As in Paper I, we have here removed one 
spatial dimension from the problem by considering vertically in- 
tegrated tori and have neglected the variations in the background 
spacetime produced by the perturbations (i.e. Cowling approxima- 
tion). Also in this case, therefore, the solution of the eigenvalue 
problem is simplified considerably and reduces to the solution of a 
single second-order ordinary differential equation. 

Overall, the results found indicate that the p-mode oscillations 
of vertically integrated tori in a Kerr spacetime share many of the 
features already encountered in a Schwarzschild spacetime. In par- 
ticular: the dependence of the fundamental eigenfrequencies on the 
position of the rest-mass density maximum and on the radial size 
of the discs, the relation of p modes with the radial epicyclic os- 
cillations of point-like particles, or the simple 2:3 harmonic se- 
quence for the lowest-order eigenfrequencies. Furthermore, this 
work extends to geometrically thick discs all of the relativistic 
disco-seismology analyses carried so far for thin discs (Okazaki 
et al., 1987; Perez et al., 1997; Silbergleit et al., 2001; Kato, 2001, 
Rodriguez et al., 2002). 

The oscillations of geometrically thick discs are important not 
only because they could produce intense gravitational radiation, but 
also because, in low-density discs, they may serve to explain the 
high-frequency quasi-periodic oscillations (HFQPOs) observed in 
low-mass X-ray binaries (LMXBs) containing a black hole candi- 
date. In these systems, in fact, the X-ray luminosity is modulated 
quasi-periodically, giving rise to distinctive peaks in the power 
spectral density which, so fare, have been found in sequences of 
small integers 2:3 (see Abramowicz and Kluzniak, 2001 for the in- 
tepretation of the results of Remillard et al., 1999 and Strohmayer, 
2001; see Remillard et al., 2002, Homan et al., 2003 for further 
observational evidences). More recently and although with a much 
poorer statistics, the possible presence of a QPO structure has been 
suggested also in the power spectrum of the light curve of the two 
brightest X-ray flares from the Galactic Center black hole (Aschen- 
bach et al., 2004). 

Using the striking analogy between the results of the numer- 
ical simulations and the observations in LMXBs, Rezzolla et al. 
(2003a) have proposed a simple model that exploits the properties 
of p-mode oscillations in thick discs and accounts for the com- 
plex phenomenology observed for HFQPOs. It is interesting to 
note while such a configuration could be produced whenever an 
intervening process modifies the Keplerian character of the flow 
near the black hole, the numerical simulations of De Villiers et al. 
(2003) have now provided a more realistic clue to how these tori 
can be generated. Similarly, the recent model proposed by Gian- 



nios & Spruit (2004) suggests a simple way in which these p-mode 
oscillations could be excited. Both of these models were not avail- 
able at the time the model by Rezzolla et al. (2003a) was proposed. 

To further develop this idea, we have here considered whether 
a purely Newtonian description of physics could be sufficient to 
account for the observations. For studying this, we have also per- 
formed a global analysis of the axisymmetric oscillations of ver- 
tically integrated non-Keplerian discs in Newtonian physics. Our 
results indicate that no major qualitative differences emerge and 
that the most important features of p modes in relativistic tori re- 
main unchanged also in the corresponding Newtonian models, al- 
beit with quantitative differences. 

The plan of the paper is as follows: in Section|2|we introduce 
the basic assumptions and equations employed in the definition of 
our general relativistic, vertically integrated tori. These equations 
will then be used to study axisymmetric oscillations both locally, in 
Section|3| and globally, in Section|4| We will first consider config- 
urations with constant distributions of specific angular momentum 
and subsequently distributions in the cylindrical radial coordinate 
distributions of specific angular momentum that increase outwards 
as power-laws. In Section [5] we comment on how these axisym- 
metric oscillations can be used to explain the HFQPOs in LMXBs 
containing a black hole and discuss the solution of the eigenvalue 
problem for vertically integrated Newtonian tori. Finally, Section|6| 
contains our conclusions. 

Hereafter, Greek indices are taken to run from to 3 and Latin 
indices from 1 to 3; unless stated differently, we will use units in 
which G — c — M© = 1. 



2 RELATIVISTIC TORI IN A KERR SPACETIME : 
ASSUMPTIONS AND EQUATIONS 

We will here assume that the torus does not contribute to the space- 
time metric, which we will take to be that external to a Kerr black 
hole. Furthermore, since we are interested in the portion of the 
spacetime in the vicinity of the equatorial plane (i.e. for values of 
the spherical angular coordinate \9 — n/2\ <C 1), we will write the 
Kerr metric in cylindrical coordinates (t,u?,<f),6) and retain the 
zeroth-order terms in the ratio (z/zu). In this case, the line element 
assumes the form (Novikov and Thorne 1973) 

ds 2 — — tU . dt 2 + -^7-(d<f> — uodi) 2 + ^—dvj 2 + dz 2 , (1) 
A uj 2 A 

where A = uj 4 + vj 2 a 2 + 2Mvja 2 , A = -ou 2 ~ 2Mm + a 2 and 
uj = 2Mazu/A. Here, M is the gravitational mass of the black 
hole and a/M is the Kerr parameter, so that the black hole angular 
momentum can be expressed as J = aM . 

The basic equations to be solved to construct models in hy- 
drostatic equilibrium are the continuity equation Va,(pu Q ) = 
and the conservation of energy-momentum, V a T a/3 = 0, where 
the symbol V refers to a covariant derivative with respect to the 
metric Q. Here, T al3 = (e + p)u a u 13 +pg al3 are the components 
of the stress-energy tensor of a perfect fluid, with u a being the 
components of the 4- velocity, p the rest-mass density, e the energy 
density and p the pressure. It is also useful to introduce an orthonor- 
mal tetrad carried by the local stationary observer and defined by 
the one-forms with components 

u = uj^/A/Adt , = VA(d(j) - udt)/w , 

oj z = dz , = VJ /\f~Adzxj . (2) 

In this frame, the components of the four velocity of the fluid are 
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denoted by it M and the 3-velocity components are defined as 



U LO a U 
«* LU%U° 



■w, z, 



(3) 



We consider the perfect fluid to follow a polytropic equation 
of state (EOS) p = kp 1 , where k and 7 = dlnp/dln p are the 
polytropic constant and the adiabatic index, respectively. Next, be- 
ing interested in a vertically-integrated description of the torus, we 
introduce a vertically integrated pressure 

P{w) = ( pdz, (4) 



and a vertically integrated rest-mass density 

cH 

E(cc) = / pdz, 



(5) 



where H = H(vj) is the local "thickness" of the torus. We further 
assume that P and E obey an "effective" polytropic EOS 



P = /CE 1 



(6) 



so that K, and r = d In P/d In E play the role of the polytropic 
constant and of the adiabatic index, respectively. As mentioned in 
Paper I, it is important to underline that {6j does not represent a 
vertically integrated polytropic EOS since the polytropic exponent 
r = r(ccr) is now a function of the position inside the torus. 

After the vertical integration, we enforce the conditions of hy- 
drostatic equilibrium and axisymmetry (i.e. assume dt = = d,/,) 
and simplify the equation of energy-momentum conservation to a 
Bernoulli-type form (Kozlowski et al., 1978) 



= — Vi m(u t ) + 



e + p 



(7) 



where I = —u^/ut is the specific angular momentum (i.e. the 
angular momentum per unit energy). Using the only relevant com- 
ponent of equations Q, we construct the equilibrium model for a 
non-selfgravitating torus in the Kerr spacetime through the force- 
balance equation 



1 



dP 



M(l-a,n) 2 /z 



■ujVl 2 



E + P dvj 



2 A/ A- A(u-Q,) 2 /z 



(8) 



where E is the vertically integrated energy density (cf. eq.|4j- 

We next perturb the hydrodynamical equations introducing 
Eulerian perturbations of the hydrodynamical variables with a har- 
monic time dependence of the type 



5V™,5V",5Q) ~ exp(-io-i) , 



(9) 



where 5Q = SP/(E + P) and where we have defined the vertically 
averaged velocity perturbations respectively as 



8V n 



1 
2H 



Sv^dz . 



2H 



— I Sv^dz . (10) 



As in our previous investigation in a Schwarzschild spacetime, we 
assume that the Eulerian perturbations in the metric functions can 
be neglected, i.e. 5g a b = (Cowling approximation; Cowling, 
1941). While this condition does not hold true in general, it repre- 
sents a very good approximation in the case of non-selfgravitating 
tori. 

Introducing now the quantities 

SU = iSV*, SW = SV^ , (11) 

to eliminate the imaginary part from the system of equations and 



after a bit of straightforward algebra, we derive the following set of 
ordinary differential equations 



A A 
o—=W + a — 

VA ^ 2 
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A 3 / 2 (Alo\ .A 3 / 2 , P' 
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sw = o. 



(12) 



sw + ln' + n 
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+ 



vj 2 A 
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Au' 
Aw\A A 



Au 

7^2 



(u - n) - 



— - 



^ A Avj a 

(13) 



gSQ + f-%SU' + \ 



p' 



E + P 



+ 



r 2 A 
A 



VJ 



Va 



6U- 



<j\TA(u - Q)Azu 2 
w^A- A 2 {lo -Q.f 



+ 



FSW = . 



(14) 



where a = l/(u t ) 2 ,f = FP/E + P, and the index ' indicates a 
radial derivative. 

Equations 1121 — 1141 are the uj- and ^-components of the per- 
turbed relativistic Euler equations and the perturbed continuity 
equation, respectively. As in Paper I, they are solved numerically 
for the eigenfrequencies and eigenfunctions of p-mode oscillations 
of an oscillating vertically integrated thick disc in a Kerr spacetime. 
Before discussing the solution of this eigenvalue problem, however, 
it is instructive to consider how wave-like perturbations propagate 
in the torus. 



3 PERTURBATIONS OF RELATIVISTIC TORI: A 
LOCAL ANALYSIS 

We now present a local analysis of the perturbed hydrodynamical 
equations 1 1 21 — i 14t which will serve to clarify the relation between 
acoustic waves and other waves that play a fundamental role in flu- 
ids orbiting in a central potential, i.e. inertial (or epicyclic) waves. 
For this analysis we assume a harmonic radial dependence in the 
perturbed fluid quantities of the type (5Q, SU, 8W) ~ exp(ifcro), 
where k is the radial wavenumber so that the local wavelength is 
A = 27r/fc. 

To make the system of equations simpler, we have removed 
from equations i 12t — i 14i those terms that are much smaller than 
the others. Details on these simplifications are discussed in the Ap- 
pendix^ The linearized perturbation equations i 12t — i 14i resulting 
from this procedure can thus be written as a homogeneous linear 
system in matrix form 



M(k,cr 




= 



(15) 
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where M is the coefficient matrix given by 
M(k,a) = 



( a A A5 

HA 

ifcf A 



n 



—2) ~\~) + r 2 (E + P) 



rVA 

A 



ifcaA\ 



(16) 



The dispersion relation is then obtained by imposing the de- 
terminant of M to be zero, thus guaranteeing that a non-trivial so- 
lution to the system of equations <15t exists. After a bit of algebra 
one obtains the dispersion relation 



2 2 . " 
a = K r H 7 

VJ 2 



zu 2 A A(lo 2 



A 



(17) 



where c s = v/ dP/ d_B is the relativistic sound speed within the 
vertically integrated disc and where the relativistic radial epicyclic 
frequency for an extended fluid object in a Kerr spacetime Ki is 
defined as 



AC 



-4 1 n- 



Alu\ 2A tn , P' 

^ + ^ (n - co) ETp 



(18) 



Here, the quantity C is just a shorthand for 



C ■: 



a(lu - n) 

w 2 A 



zu 2 A 
A 



Au' 



Aiv 
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+ 



A 2 loco' QA 2 lj' 
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, _ ZD ( AUJ 

A \zu 2 



(19) 

Equations J 171 and J 1 81 include corrections coming from the 
rotation of the black hole and pressure gradients from the fluid dis- 
tribution and reduce respectively to the dispersion relation and to 
the relativistic epicyclic frequency for an extended object obtained 
in Paper I for a Schwarzschild black hole (cf. eqs. 42 and 43 of Pa- 
per I). When the motion is almost Keplerian, the pressure gradients 
are negligible and expression i 1 8i reduces to the expression for the 
radial epicyclic frequency for point-like particles (Okazaki et al., 
1987) 



'Kcp 



M_ 

777 



- 6M/zu ± 8a^M/vu 3 - 3a 2 /vu 2 
(1 ± Sa^M/vu'*) 2 



(20) 

where the ± sign distinguishes orbits that are corotating from those 
that counter-rotating relative to the black hole spin. 

The behaviour of the radial epicyclic frequency 1201 is shown 
in Fig.Qwhere different curves refer to different black hole spins 
and where we consider corotating orbits only. For each curve, fur- 
thermore, a dashed line refers to the "increasing" branch of the 
epicyclic frequency {i.e. with kJ. > 0), while a solid line refers to 
the "decreasing" branch (i.e. with k' v < 0). The solid circles indi- 
cate the locations of the maximum frequencies for each black hole 
spin and will be used subsequently in Fig.|3| Finally, the dotted line 
indicates the Newtonian radial epicyclic frequency which will be 
used for the discussion in Section [STTI The numerical values refer 
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Figure 1. Radial epicyclic frequency for different values of the black hole 
spin. In each curve, a dashed line refers to the "increasing" branch (i.e. with 
kJ > 0), while a solid line to the "decreasing" branch (i.e. with k[ < 0). In- 
dicated with solid circles are the locations of the maximum frequencies for 
each black hole spin, while the dotted line refers to the Newtonian values 
for the radial epicyclic frequency. The numerical values have been com- 
puted for a black hole with mass M = IOMq and the radial extents have 
been expressed in units of gravitational radii r g = GM/c 2 . 



to a black hole with mass M = IOMq and the radial extents have 
been expressed in units of gravitational radii r g = GM/c 2 . 

Overall, the dispersion relation J 171 shows that the propaga- 
tion of small perturbations in a fluid rotating around a Kerr black 
hole is characterized by two main features. The first one is of purely 
acoustic nature (i.e. with a oc kc a ); the second of is of purely iner- 
tial nature (i.e. with u a k,) and is reminiscent of the oscillations 
that a particle orbiting in a central potential experiences when per- 
turbed. For these reasons, the resulting waves are usually referred 
to as inertial-acoustic waves and the corresponding modes can be 
classified as p modes. The determination of the eigenfrequencies 
and eigenfunctions of these modes will be discussed in the follow- 
ing Section. 



4 PERTURBATIONS OF RELATIVISTIC TORI: A 
GLOBAL ANALYSIS 

We now turn to a global analysis of the axisymmetric oscillation of 
relativistic tori in a Kerr spacetime by solving the system of equa- 
tions CTH-fHl as an eigenvalue problem. In particular, the solution 
is found using a "shooting" method (Press et al., 1986) in which, 
once the appropriate boundary conditions are provided, two trial 
solutions are found, starting from inner and outer edges of the disc 
respectively, and these are then matched at an intermediate point 
where the Wronskian of the two solutions is evaluated (see Paper I 
for details). This procedure is iterated until a zero of the Wronskian 
is found, thus providing a value for a and a solution for 5Q, SU, 
and SW. The numerical methods employed here to solve the eigen- 
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Figure 2. Schematic representation of the position of the relevant radii in 
the cases of "constant" and "power-law" distributions of specific angular 
momentum. In particular, while ro C usp and ro max indicate the positions of 
the cusp and of the rest-mass density maximum, ro m b and ro ms indicate 
the positions of the marginally bound and marginally stable orbits. 
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Figure 3. Possible locations of the rest-mass density maximum ro max for 
tori with constant specific angular momentum distribution and for all possi- 
ble black hole spins. The solid lines refer to the smallest and largest possi- 
ble positions of ra max , while the dashed line indicates the locations of the 
maximum of epicyclic frequency k t (cf. the solid circles in Fig. HI. 



value problem are the same as those discussed in Paper I, where a 
more detailed discussion can be found. 

It is useful to briefly summarize how the background model 
for the torus is constructed (see also Paper I) as this will help clar- 
ifying the discussion of the results presented in the following Sec- 
tions. First of all, a value of the black hole spin a/M is chosen 
and the distribution of the specific angular momentum I = t(zu) 
is selected. The positions of the cusp and of the maximum rest- 
mass density in the torus are then obtained by imposing that the 
specific angular momentum at these two points coincides with 
the Keplerian value (see Fig. [2] for a schematic view). The inner 
edge of the torus zui n is determined by fixing the potential gap, 
AW in = Win - Wcus P , defined as 

AW in = ln[(-ut) fa ] - ln[(-« t )cusp] - / , (21) 

with AWi n = corresponding to a torus filling its outermost 
equipotential surface. Next, the hydrostatic balance equation JSJ is 
integrated from -coin to the outer edge of the torus tu ou t, defined 
as the position at which P = 0. Sequences of tori having different 
radial extents for a given distribution of specific angular momen- 
tum and black hole spin can then be constructed by varying the 
potential gap A Win. In such sequences, the tori will have the same 
rest-density maxima n7 max . 

Clearly, different distributions of specific angular momentum 
will produce tori with different positions of the cusp and of the 
maximum rest-mass density. In order to investigate how the ax- 
isymmetric oscillations depend on these properties, we have con- 
structed models with different distributions of specific angular mo- 
mentum, taken to be either constant within the torus, or increasing 
outwards according to a generic power-law in radius (cf. eq. l24> . 
The main properties of the various models considered are summa- 
rized in Table Q where they are presented in terms of dimension- 



less quantities. A simple conversion of the dimensionless eigenfre- 
quency a and of the specific angular momentum I is obtained, for 
instance, through the following relations 

a = a {wj {—) c ' (22) 

and 

4.1 Constant specific angular momentum tori 

4.1.1 Fixed black hole spin 

We first consider tori with a distribution of specific angular mo- 
mentum that is constant in space, i.e. l(ru) — const., not only 
because this choice is mathematically simpler to investigate, but 
also because it allows for a straightforward determination of the 
permitted range of equilibrium models. In this case, in fact, for 
a torus of finite size to exist, the value of the specific angular 
momentum must satisfy ^ ms < I < l m t, where ^ ms and l! m b are 
the specific angular momenta of the marginally stable and of the 
marginally bound orbit, respectively. For a Schwarzschild black 
hole, £ ms = 3\/6/2 ~ 3.67 and ^ m b = 4, while they depend of the 
the black hole spin for a Kerr black hole. Because the position of 
the rest-mass density maximum iu max will depend on the specific 
value of I chosen, it is then possible to define the range in which 
Umax lies. In other words, all of the possible finite-size tori will 

have (tCTmax) 

MIN ^ ^max < (^maxJmAX ' ^ defined 

(^max) MIN as the value of cf7 max when i = £ ms , and (cf7 m ax) MAX 
as the value of G7 max when I = £ mb . Clearly, (n7 max ) MIN also 
coincides with the location of the marginally stable orbit, while 
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Table 1. Main properties of the equilibrium models studied. From left to right the columns show: the black hole spin parameter.the type of specific angular 
momentum distribution, the constant coefficient l c and the power-law index q (cf. eg. 1241 . the position of the maximum density point m max , the position of 
the inner and outer radii of the torus ro in and zu out , the position of the maximum density point with respect to the epicyclic frequency curve, and the potential 
gap at the inner edge of the torus AVV; n . This latter quantity is reported because not all equilibrium models reported in this table have been constructed as the 
ones filling the largest closed equipotential surface (i.e. P16-P19). 
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(■ci7max) MA x provides the position of the rest-mass density maxi- 
mum for a torus whose outermost equipotential surface closes at 
infinity. 

Furthermore, because in a Kerr spacetime the values of £ms 
and £ m b will depend on the black hole spin a/M, it is possible to 
define the set of all the possible values that tu max can assume for 
constant angular momentum tori construncted in a Kerr spacetime. 
This is shown as a shaded area in Fig.|5| where we have also plotted 
with solid lines the different values of (ro ma x) MIN and (sj m a X ) MAX 
as the black hole spin is increased from to 1 1 . Also shown with 
a dashed line are the locations of the maximum of epicyclic fre- 
quency for different black hole spins (cf. the solid circles in Fig.0. 

Note that for tori with constant specific angular momentum 
distributions around rapidly rotating black holes with spin parame- 
ter a/M larger than a*/M (cf. Fig.[5), all of the possible values of 
Umax will be located on the "increasing branch" of the epicyclic 
frequency (i.e. at positions where k' t > 0). As a result, the closer 
the maximum density point is to the black hole horizon, the smaller 
the value of the epicyclic frequency at that point (cf. eq. l20> . For ro- 
tations of the central black hole smaller than a,/M, on the other 



1 We have not considered here counterotating black holes, i.e. with a < 0. 



hand, the models with tu max in the "decreasing branch" of the 
epicyclic frequency (i.e. at positions where k! x < 0) are also al- 
lowed together with those on the increasing branch. In this case, if 
the maximum rest-mass density point coincides with the location 
of the maximum point of the epicyclic frequency (i.e. at positions 
where k[ = 0), the eigenfrequency of the fundamental mode for a 
given radial extent will have the largest value for any a < a*. The 
locus of these points is indicated with a dashed line in Fig.|3| 

In order to investigate the dependence of the fundamental 
mode of oscillation on the position of the maximum rest-mass den- 
sity, we have first considered sequences of tori having different ra- 
dial extents L = vj out — vj- ln . All sequences are built for the same 
value of the black hole spin, while each sequence is characterized 
by a specific value of ro maj( , which is chosen to be either on the in- 
creasing or on the decreasing branch of the corresponding epicyclic 
frequency. The results for a black hole spin a/M — 0.5 are sum- 
marized in the two panels of Fig.|4] which show the values of the 
eigenfrequency of the fundamental mode of oscillation / as a func- 
tion of the disc radial extents L. 

In the left panel, in particular, we show the results for mod- 
els whose maximum rest-mass density is located on the decreasing 
branch of the radial epicyclic frequency. As expected for modes 
behaving effectively as sound waves trapped in a rotating fluid, the 
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Figure 4. Eigenfrequencies for the fundamental mode of axisymmetric p modes for tori with a constant distribution of specific angular momentum. Each line 
corresponds to a sequence of tori having the same ro max but different radial extents L and the solid circles correspond to the values of the Keplerian radial 
epicyclic frequency 1201 at ro max . Both panels refer to a Kerr black hole with a/M = 0.5, but while the right panel has ro max located in the increasing 
branch of the epicyclic frequency (i.e. rej. > 0), the left one has ro max located in the decreasing branch of the epicyclic frequency (i.e. rej. < 0). 
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Figure 5. Eigenfrequencies for the fundamental mode of axisymmetric p modes for tori with a constant distribution of specific angular momentum. Each 
line refers to a sequence of tori for the same black hole spin a / M, but having different radial extents L and the solid circles correspond to the values of the 
Keplerian radial epicyclic frequency 1201 at ro max . In each panel, the sequences have been computed keeping fixed the position of ra max , which are however 
different according to whether the black holes are slowly rotating (left panel) or rapidly rotating (right panel). 



eigenfrequencies decrease like L _1 as the radial extent of the torus 
increases. Note also that, as was shown in Paper I for tori orbit- 
ing around a Schwarzschild black hole, the eigenfrequencies of the 
fundamental mode tend to the values of the radial epicyclic fre- 
quency at t*7 max as the radial dimension of the discs tends to zero 



(In Figs.|4]0and[9]the filled circles show the values of the Keple- 
rian radial epicyclic frequency 1201 at ro max .). These results confirm 
the expectation that, as their size diminishes, the discs effectively 
behave as rings of particles in circular orbits and have as their char- 



8 Montero, Rezzolla and Yoshida 



acteristic frequency the epicyclic frequency at the maximum rest- 
mass density point. 

Note also in the left panel of Fig.[4]fhat models with the same 
radial extent and different constant angular momentum distribu- 
tions, i.e., tori with different locations of ?i7 max , have fundamental 
eigenfrequencies that become larger as the rest-mass density max- 
imum approaches the black hole. This behaviour can be explained 
simply since, for L — *• 0, the fundamental eigenfrequencies con- 
verge to the values of the epicyclic frequencies and the latter in- 
crease for smaller radii as one moves along the decreasing branch. 

The opposite behavior is observed in the right panel of Fig. [4] 
but in this case the maximum density point of all models is located 
on the increasing branch of the epicyclic frequency. As a result, 
models with the same radial extent have fundamental frequencies 
that become smaller as the rest-mass density maximum approaches 
the black hole. Therefore, for any given radial extent, the model 
with the largest fundamental-mode eigenfrequency will be the one 
with G7 max located at the maximum of the epicyclic frequency (i.e. 
along the dashed line in Fig.[3j- 

4.7.2 Influence of the black hole spin 

In order to investigate the influence of the spin of the central black 
hole on the axisymmetric oscillations, we have solved the eigen- 
value problem for sequences of tori having the same location of 
ro max but different black hole rotation rates. 

We note that, in the case of constant specific angular momen- 
tum, it is not possible to construct sequences of tori having the same 
location of the rest-mass density maximum for all the possible val- 
ues of the black hole spin. This is due to the fact that for any given 
ro m „, the existence of marginally stable and marginally bound or- 
bit limits the range of values that a/M can assume. This is quite 
apparent in Fig.|3|where, for any chosen value of m max , the black 
hole spin cannot span the whole range from to 1. 

Because of this constraint, we show in the left panel of Fig. [5] 
sequences of tori having the same value of tu max /A/ = 7.0 and 
black hole spins varying from a/M — 0.0 to a/M — 0.5. Clearly, 
these sequences share many of the properties discussed for Fig.|4| 
and, in addition, for any given radial size of the torus, the frequency 
of the fundamental mode will be larger for larger black hole spin. 
Indeed, this is what is expected, since the radial epicyclic frequency 
for a point-like particle (tr) Kep increases monotonically with a for 
any value of va (cf. Fig.0. The same behaviour is observed also 
when tu max /M = 4.5, which allows for larger rotations rates of 
the black hole and is shown in the right panel of Fig. [5] 

4.1.3 Harmonic Sequence 

An interesting and potentially important feature of axisymmetric 
p-mode oscillations of tori with constant specific angular momen- 
tum is that the lowest-order eigenfrequencies appear in a sequence 
2:3:4..., with a good accuracy (i.e. within 10%) and rather indepen- 
dently of the details of the disc properties. This feature was first dis- 
covered through the fully nonlinear numerical simulations of Zan- 
otti et al. (2003) and subsequently confirmed through a perturbative 
analysis in a Schwarzschild spacetime in Paper I, and through ad- 
ditional independent numerical simulations (Kluzniak et al., 2004; 
Lee et al., 2004). 

A harmonic relation for the lowest-order eigenfrequencies 
continues to exist also for tori with constant specific angular mo- 
mentum in a Kerr spacetime and rather insensitively of the black 
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Figure 6. Ratio of the first overtone and the fundamental mode of os- 
cillation for tori with constant specific angular momentum distributions. 
The different sequences refer to the different models described in Table HI 
with the solid lines referring to ro max = 7.0 r g and the dashed ones to 
TOmux = 4.5 r g . 

hole spin. Fig.|6|shows the ratio between the frequencies of the first 
overtone oi and the fundamental mode / for some representative 
tori with constant specific angular momentum distributions consid- 
ered in Table Q Different line types correspond to sequences that 
have been constructed with different locations of the rest-mass den- 
sity maximum. All of the sequences span values of a/M from to 
0.8 without showing major qualitative differences. 

Since it is not the result of a mathematical constraint but, 
rather, the consequence of a global mode of oscillation, the 2:3 ra- 
tio between the fundamental mode and its first overtone is satisfied 
with an accuracy that is in general within 10%, but is, nevertheless, 
not exact. A number of different elements can contribute to a devi- 
ation from an exact relation among integers and the results of our 
calculations show that conditions such as the size of the disc, the 
location of the rest-mass density maximum, the black hole spin, the 
distribution of specific angular momentum, or the EOS considered, 
can all influence this departure. Note also that in Fig.|6|some of the 
models (e.g. C4b and C3a) appear to have a rather large variation in 
the ratio oi//. It should be recalled, however, that these tori exist 
only over a very small range in the radial extent L (either because 
of the high value of the specific angular momentum or because of 
the large spin of the black hole). Being trapped in such small tori 
rather close to the black hole horizon where the angular velocity 
is larger, the p-mode oscillations will be particularly sensitive on 
variations of the radial size. 

Another interesting feature of the lowest-order eigenfrequen- 
cies is that, for any given black hole rotation speed, their ratio tends 
to a single value oi // ~ 1.52 as the size of the tori tends to zero (cf. 
Fig. [6). This is in good qualitative agreement with what was found 
in the perturbative analysis of slender tori (Blaes 1985), where the 
dispersion relation for generic perturbations of Newtonian slender 
tori having constant distributions of specific angular momentum, 
was derived (cf. eq. 1.8 of Blaes 1985, with j = 1, and |fc| = 0). 
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4.2 Non-constant specific angular momentum tori 

We next investigate axisymmetric p-modes in tori with non- 
constant distributions of specific angular momentum. Clearly, in 
this case the problem is highly degenerate since there is a complete 
arbitrariness on what is the most appropriate distribution of angu- 
lar momentum. The first guide in the infinite possible choices is 
that they should correspond to stable fluid configurations. While in 
Newtonian theory one can make use of the simple Rayleigh stabil- 
ity criterion for rotating inviscid fluids dl/dzo > (Tassoul 1978), 
the situation in a general relativistic framework is more complex. 
In particular, the condition of stability requires that, for the perfect 
fluids considered here, the gradient of the specific angular momen- 
tum must never point toward the interior of the quasi-circular level 
surfaces of I (Seguin 1975, Abramowicz and Prasanna 1985). Be- 
cause of this, we have used a simple prescription for the specific 
angular momentum in which it increases outwards with a power- 
law distribution of the type 

1 = Ic-uj" , (24) 

where both £ c and q are positive constants (cf. Table Q for the dif- 
ferent values used). 

Following the procedure presented in the previous Section, we 
have constructed tori by selecting the angular momentum distribu- 
tion such that the position of the cusp is always located between 
the marginally bound and the marginally stable orbits. Note that, in 
principle, this is not the only possible choice. Non-constant distri- 
butions of specific angular momentum, in fact, allow for the con- 
struction of tori having the cusp at radii larger than the position 
of the marginally stable orbit. In such tori, matter lost through the 
cusp as a result of small perturbations would be able to find new 
stable circular orbits and, as a result, the accretion onto the black 
hole will require mechanisms other than simple perturbations in the 
flow. Here, we have not investigated the properties of these discs 
whose dynamics, however, has been considered in detail by Zanotti 
et al., (2004). 

Overall, the solutions for tori with non-constant distributions 
of specific angular momentum share qualitative similarities with 
those found when £ = const. In particular, the main properties of 
the p-mode oscillations in this case can be summarized as follows: 
Firstly, the eigenfrequencies depend on the location of the rest- 
mass density maximum K7 max and on whether the latter is located 
on the increasing or decreasing branch of the radial epicyclic fre- 
quency. Secondly, all eigenfrequencies tend to the radial epicyclic 
frequency at the location of ii7 max in the limit of small tori. Thirdly, 
as the result of the presence of an evanescent-wave region (see Pa- 
per I for a detailed discussion), the eigenfunctions become vanish- 
ingly small at the inner edge of the disc, while they do not do so at 
the outer edge. Finally, for any sequence built with fixed ro max , the 
eigenfrequencies increase as the spin of the black hole is increased. 

Despite these similarities, a difference does emerge when in- 
vestigating non-constant angular momentum tori. In this case, in 
fact, the ratio between the first overtone and the fundamental fre- 
quency of oscillation Oi/f has a behaviour which is more complex 
than the one found when t = const. In particular, it depends not 
only to the position of ro max , but also on the values of i c and q, as 
well as on a/M. Most importantly, while the eigenfrequencies re- 
main in a harmonic sequence for sufficiently large discs, departures 
from a 2:3 ratio can be observed for very small discs. This is shown 
in FigQ where the oi/f ratio is plotted for sequences of models 
orbiting around central black holes rotating at different rates and 
having different distributions of specific angular momentum (see 



o 



1 


5 


ii 


i i i i i i i i 

l = power — law _ 
a = 9.0 

max 


1 


4 


P10 


P14 P17 - 

^ — P19 _ 




1 


3 


-"i i i i lii 


i i i i i lii 



10 20 30 



Figure 7. Ratio of the first overtone and the fundamental mode of oscillation 
for tori with a power-law distribution of specific angular momentum. The 
different sequences distinguish different choices of the index q, but all refer 
to the same value for £ c (cf. TablefTl. 

Table Q for the model descriptions). In all the cases considered, 
however, the departure from a 2:3 ratio is less than 20%. 

Finally, note that, as the size of the discs tends to zero, the 
ratios do not tend to the same value, as was the case for con- 
stant specific angular momentum tori but, rather, to a value that 
depends on the index q of the power-law. This is shown in Fig. [H] 
where we illustrate results for sequences of tori with the same l c 
but different indices q, built around the same Kerr black hole with 
a/M = 0.5. Also in this case, a closer agreement with a harmonic 
relation between the lowest-order eigenfrequencies is recovered for 
sufficiently large discs. 



5 IMPLICATIONS FOR HFQPOS IN LMXBS 

Observations of X-ray emissions from binary systems comprising 
a stellar-mass compact object have long been considered important 
tools to test General Relativity in strong-field regimes. In particular, 
the high frequency quasi-periodic oscillations (HFQPOs) observed 
in LMXBs binaries containing a black hole candidate have been 
proposed as a tool for measuring black hole properties such as mass 
and spin in a direct way. In such systems, the X-ray observations 
show a particularly rich phenomenology in which the luminosity 
is modulated quasi-periodically, giving rise to distinctive peaks in 
the power spectral density. A puzzling feature of the QPO frequen- 
cies in these systems is that they are found in sequences of small 
integers (i.e. 1:2, 2:3, or 1:2:3; Abramowicz and Kluzniak, 2001; 
Remillard et al., 2002; Homan et al., 2003). 

Numerous models have been suggested to explain the HFQ- 
POs and their peculiar harmonic property (see Abramowicz and 
Kluzniak, 2004 for a recent review) but the scarcity of the obser- 
vational data, on the one hand, and the crudeness of the different 
models, on the other hand, have not yet made it possible to deter- 
mine which model provides the most likely description of the pro- 
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Figure 8. Ratio of the first overtone to the fundamental mode shown as 
a function of the radial extent of the torus L. The tori are modelled with 
a power-law distribution of specific angular momentum and the different 
lines refer to different choices of the index q, but all assume the same value 
of£ c = 2.73. 



cesses responsible for the QPOs. It should be noted the first steps 
towards a more quantitative and physically realistic description of 
the emission properties of possible QPO models have recently been 
taken (Schittman and Bertschinger, 2004; Schnittman, 2004). 

Among the proposed models, one is particularly simple as it 
is based on the assumption that the accretion disc around the black 
hole terminates with a sub-Keplerian part, i.e a torus of small size 
(Rezzolla et al., 2003a). This model was motivated by the numer- 
ical simulations of relativistic tori in a Schwarzschild spacetime 
(Zanotti et al., 2003) and by the subsequent perturbative study of p- 
mode oscillations of relativistic vertically integrated tori presented 
in Paper I. A Cornerstone of this model is the evidence, both nu- 
merical and analytic, that in these objects the fundamental mode 
and the first overtones are found to be in the harmonic sequence 
2:3:4... to a good precision and in a very wide parameter space. 

Because the p modes considered in the above model represent 
the basic oscillation modes of a non-Keplerian disc, some of their 
features are not expected to be very sensitive to the properties of 
the spacetime which determine the motion of the fluid. The har- 
monic 2:3 relation between the fundamental frequency and the first 
overtone is one such feature and, indeed, this property has been en- 
countered essentially unmodified in Schwarzschild and Kerr space- 
times, but also in a Newtonian gravitational potential. Because of 
this, it is then natural to consider whether a fully Newtonian de- 
scription of the physics is sufficient to account for the phenomenol- 
ogy observed in HFQPOs, or whether General Relativity is indeed 
necessary. This question in addressed in the following Section. 

5.1 Can Newtonian Physics explain the HFQPOs? 

In order to assess whether p-mode oscillations in purely Newto- 
nian tori could account for the HFQPO phenomenology, we have 
performed a global analysis of the axisymmetric oscillations of ver- 



tically integrated non-Keplerian discs in Newtonian physics and 
compared the results with those obtained in a general relativistic 
framework. 

Consider therefore a non-selfgravitating perfect fluid torus in 
equilibrium and orbiting a central object of mass Mo producing a 
gravitational potential ^1/. As in the general relativistic case, the sys- 
tem of equations describing the torus will be simplified by consider- 
ing vertically integrated or vertically averaged quantities, which we 
indicate as P and E, in analogy with the corresponding relativistic 
quantities (cf. eqs.|4]and|5j- In order to handle simpler expressions 
we normalize all quantities as 
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where i — w, cj>, while Eo is the maximum mass density and where 
we have defined 
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(26) 

The equilibrium configurations are then obtained by solving 
the hydrostatic equilibrium equation which, in terms of the normal- 
ized quantities, is written simply as 



E 



(27) 



where ^ = $?/Co and SI = fi/fio, with fio being the angular 
velocity at zuq. Introducing now the Emden function B and its re- 
lated quantity Q = NO, the normalized rest-mass density E can 
written as E = O^, and the equation of hydrostatic equilibrium is 
expressed as an ordinary differential equation for Q, i.e. 



Q' — Q — . 



(28) 



Once the angular velocity distribution Q = f2(w) is fixed and a 
choice is made for the inner radius zu- ln , equation <28> can be inte- 
grated to provide the background Newtonian torus in equilibrium. 

For simplicity we now restrict our attention to constant spe- 
cific angular momentum distributions so that the equation of hy- 
drostatic equilibrium I28i simplifies to 



/ e 



(29) 



Before discussing the solution of the eigenvalue problem, it is use- 
ful to recall the features of the Newtonian models that represent 
important differences from the relativistic counterparts. The first 
one is that the right-hand side of equation I29i is just the derivative 
of an effective potential (i.e. the one in the round brackets) and it is 
therefore straightforward to calculate the inner edge of the largest 
possible torus as the radial position at which this potential vanishes, 
i.e. (ro in ) MIN = £ 2 /2. 

The second feature to notice is again easily derived from equa- 
tion <29t and shows that the position of the mass density maximum 
(i.e. Q' = 0) is simply given by n7 max = £ 2 , coinciding with the 
position at which the specific angular momentum assumes a Keple- 
rian value. As in relativistic tori, S mal grows with increasing spe- 
cific angular momentum; unlike relativistic tori, however, no lower 
or upper limits exist for 5j ma x- Indeed, for vanishingly small values 
of t, it is possible to construct tori whose pressure maxima lie arbi- 
trarily close to the rotation axis. This property is the consequence 
of the fact that, in Newtonian physics and for constant distributions 
of specific angular momentum, I is not constrained to lie in a finite 
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Newtonian tori with constant angular momentum distribution orbiting a 
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interval. This represents a major difference with respect to the gen- 
eral relativistic picture and, as will be discussed later on, makes the 
parameter space for the background models effectively unbounded. 

We next proceed to consider the eigenvalue problem and, as 
we did for the general relativistic case, we neglect perturbations in 
the gravitational potential (Cowling approximation) and introduce 
perturbations in the velocity and pressure with a harmonic time 
dependence of the type ~ exp(— iat). To simplify the notation, 
hereafter we will omit the "tilde" on all variables, which however 
should be meant to refer to normalized quantities. As a result, the 
perturbed Euler and continuity equations can be written as 
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(32) 



where, in analogy with the corresponding relativistic quantities, SU 
and SW are the vertically averaged radial and azimuthal velocities, 
respectively (cf. eqs. HOI and ll 1> . 

After some simple algebra, equations J30I — J32I are cast into a 
single, second-order ordinary differential equation 
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for the quantity SQ defined as 
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SQ' 
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ND 



5Q = 0, (33) 



(34) 



- 2ft(2ft + vjQ!) 



(35) 



measures the deviation of the eigenfrequency a from the Newto- 
nian radial epicyclic frequency («?) Newt = 2Q(2Q + zjQ'), which 
coincides with the orbital frequency when the orbital motion is Ke- 
plerian, i.e. (/t r ) Nowt = Q Kc . Note that while equation 1331 ap- 
pears to be singular at the radial position where D = 0, it is actually 
regular everywhere since SQ' and D vanish at the same position (cf. 
eas.l3llandl32l 

The eigenfrequencies and eigenfunctions of the axisymmet- 
ric p modes described by equations <30t - <32| have been computed 
using the same numerical method discussed in the previous Sec- 
tions and some representative results are presented in Fig. [5] More 
specifically, we show the eigenfrequencies of the fundamental 
mode of oscillation for three sequences of constant angular momen- 
tum Newtonian tori orbiting a central source of gravitational poten- 
tial with a mass of IOMq . The sequences have been built keeping 
the position of the density maximum fixed (i.e. £7 max = 4.5, 7.0 
and 9.0, respectively) and varying the position of the inner radius to 
produce tori of different radial sizes. The inset, on the other hand, 
shows the ratio of the first overtone to the fundamental eigenfre- 
quency as a function of the radial size of the torus for the three 
different sequences. A first glance at Fig.|5|and a comparison with 
Figs.|4||5|and|6| shows that no major qualitative differences emerge 
in a Newtonian framework and that all of the most important fea- 
tures of p modes in relativistic tori remain unchanged also in the 
corresponding Newtonian models. Most notably, the fundamental 
eigenfrequencies depend on the position of the mass density and 
on the radial sizes of the discs, increasing as the latter decrease. In 
addition, as the radial sizes tend to zero, the fundamental eigenfre- 
quencies tend to the values of the radial epicyclic frequencies at the 
position of the mass density maxima. Finally, along all sequences, 
the first and second eigenfrequencies are in a harmonic sequence 
2:3 of small integers. 

Besides this qualitative agreement, however, quantitative dif- 
ferences are also present and can be easily appreciated by com- 
paring, for instance, the curve for a/M — in the left panel of 
Fig.|5|with the dotted line in Fig. [9] Indeed, the simultaneous pres- 
ence of qualitative similarities and quantitative differences between 
Newtonian and general relativistic models is not surprising and re- 
flects the fact that the p modes are fundamental modes of oscillation 
of orbiting fluid objects. As such, they depend only in part on the 
gravitational field (or background spacetime) in which the fluid is 
moving. Indeed, also for rotating stars, the qualitative features of 
the p-mode oscillations are preserved when going from a Newto- 
nian description to a full general relativistic one, with quantitative 
differences being present, however, and depending on the strength 
of the spacetime curvature. 

Finally, note that because the radial epicyclic frequencies rep- 
resent the asymptotic limit of the fundamental eigenfrequencies 
for vanishing torus sizes and since the Newtonian epicyclic fre- 
quency is always larger than the corresponding general relativis- 
tic one for all values of the black hole spin (cf. solid and dot- 
ted lines in Fig. 0, all of the Newtonian frequencies along a 
sequence with £ = const, will be larger than the correspond- 
ing relativistic ones. This, together with the fact that the Newto- 
nian epicyclic frequency diverges for vanishingly small radii [i.e. 
( K r) N(!Wt = ^Kop ~ 3y/2 ], makes it clear that for any observed 
QPO frequency it will be possible to find a Newtonian torus of size 
L that could produce p mode oscillations with the correct frequen- 
cies and in a 2:3 ratio. Stated differently, while the parameter space 
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in the (/, L) plane for relativistic tori is bounded above by the se- 
quence approaching the largest epicyclic frequency admissible by 
a Kerr black hole with spin a/M — 1, this restriction does not ex- 
ist for Newtonian tori and a Newtonian physicist, who ignores the 
existence of marginally stable orbits or of an event horizon, would 
be able to find an equilibrium solution at every point in that plane. 

As a result, and at least in principle, Newtonian physics could 
explain the HFQPOs if the only information available is the mass of 
the central object M and the frequencies of the modulation in the 
X-ray luminosity a. The only way to remove this ambiguity and 
thus disprove the interpretation of the Newtonian physicist, would 
be to make combined measurements of M and a, together with the 
radial extension of the torus L, or of the position of the mass density 
maximum cf7 max . While both these observations wwould probably 
be difficult to make in practice, similar difficulties are found also in 
proving the existence of stellar-mass black holes in X-ray binaries 
(Abramowicz, Kluzniak and Lasota, 2002). This is simply the con- 
sequence of the fact that no limit exists in Newtonian physics to the 
compactness of the sources of gravitational potential. 



6 CONCLUSIONS 

We have presented the study of the oscillation properties of non- 
selfgravitating tori orbiting around Kerr black holes. This work ex- 
tends the previous investigations of the oscillation properties of rel- 
ativistic tori in a Schwarzschild background and their importance in 
explaining the HFQPOs in X-ray binaries containing a black hole 
candidate (Rezzolla et al., 2003a). Following the same approach 
presented in these previous works, we have here considered the ax- 
isymmetric p— mode oscillations of relativistic tori assuming a ver- 
tically integrated description to reduce the eigenvalue problem to 
the solution of a system of coupled ordinary differential equations. 

We have first considered a local analysis in a Kerr spacetime 
and determined the relations between acoustic and inertial waves, 
showing that both are present in the oscillations of geometrically 
thick discs and that they play a different role depending on the 
radial size of the disc and on the position of the rest-mass den- 
sity maximum. We have then computed the eigenfunctions and the 
eigenfrequencies of the axisymmetric p modes for a large variety 
of background models which differed either in the spin of the black 
hole or in the distributions of specific angular momentum consid- 
ered. The latter, in particular, have been considered to be constant 
or to increase outwards following a power-law in radius. 

On the whole, the p-mode oscillations of vertically integrated 
tori in a Kerr spacetime possess all of the most important features 
already encountered in a Schwarzschild spacetime. Firstly, the fun- 
damental eigenfrequencies depend on the position of the rest-mass 
density maximum and on the radial size of the discs, increasing 
as the latter decreases. Secondly, the fundamental eigenfrequencies 
tend to the values of the radial epicyclic frequencies at the posi- 
tion of the mass density maxima as the radial sizes of the tori tend 
to zero. Finally, for tori constructed with constant distributions of 
specific angular momentum, the first and second eigenfrequencies 
are in a harmonic sequence 2:3 of small integers. This ratio is not 
exact but is very accurately satisfied with deviations of ~ 10% at 
most. For non-constant distributions of angular momentum, on the 
other hand, the harmonic sequence is still present for sufficiently 
large tori, but the deviation from a precise 2:3 ratio increases, be- 
ing ~ 20% at most. 

We have also investigated the implications that p-mode os- 
cillations could have on the HFQPOs observed in X-ray binaries 



containing a black hole candidate. In a model recently proposed, 
the presence of a sub-Keplerian flow (i.e. a torus) near the black 
hole and its oscillations are used to explain many of the features 
of the HFQPOs and, in particular, the 2:3 ratio in the peaks of the 
power spectral density (Rezzolla et al., 2003a). In this model, the 
measure of the oscillation frequencies and the knowledge of the 
p-mode properties can be used to measure directly the black hole 
properties, such as the mass and spin, once a reasonable estimate 
of the torus size is made. 

In view of this, we have considered whether an equivalent in- 
terpretation of the HFQPOs phenomenology can also be made in 
terms of a purely Newtonian description of physics. To this scope, 
we have performed a global analysis of vertically integrated Newto- 
nian tori following the same mathematical and numerical approach 
developed for tori in a Kerr spacetime. Our results indicate that, at 
least qualitatively, the p-mode oscillations of Newtonian tori pre- 
serve all of the properties encountered in their relativistic counter- 
parts. Furthermore, while quantitative differences do exist between 
the two frameworks, it may be difficult to distinguish between the 
two unless the oscillation frequencies are measured together with 
the radial extension of the torus L, or through the measurement of 
the position of the mass density maximum o7 max . 
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first and the second one, for which the radial derivative introduces 
a factor ik balancing the smallness of F. 



APPENDIX A: ON THE DISPERSION RELATION FOR 
TORI AROUND KERR BLACK HOLES 

We here briefly discuss the simplifications of equations J12> - J14t 
that lead to the homogeneous linear system with matrix form 1151 . 
In particular, we start from the analogy of equation 1121 with 
its Newtonian counterpart, equation <3U . whose pressure-gradient 
term SQ' and centrifugal term 2-cuflSW are of about the same or- 
der. The reason for this is that the inertial-acoustic modes are the 
results of small deviations from the equilibrium between the cen- 
trifugal force and the pressure gradients, which produce the small 
perturbative radial velocity field expressed by the first term in equa- 
tion OTi . 

We expect this to be true also in the general relativistic case 
and assume, therefore, that the second and the third terms in equa- 
tion J 121 are of the same order. We also note that since it* = 0(1), 
then a — 0(1) and thus the leading order coefficient in the square 
bracket multiplying SW is the first one, i.e. (A 3/2 / A)(A/w 2 )'Q.. 

Because now A ~ zu 2 and A ~ zu 4 , it is not difficult to show 
that the ratio of the third to the first term in equation 1 1 21 is given 
by 

a{A/zu 2 )SQ' k SQ 

T777 j (Al) 



(A 3 / 2 /A)(A/zu 2 )'QSW nSW' 

where we have replaced radial derivative with a factor ik. Since we 
assume the ratio IA1I to be 0(1), we can deduce that 

Using now this order of magnitude estimate, we can consider 
the azimuthal Euler equation d 1 3 i and evaluate the ratio between 
the first and third terms, which is given by 

(azu 2 ^A/A)SW 1 SW k 

(A3) 



[Aaa(uj- fl)/Azu 2 ]8Q fi SQ Q 2 ' 

As a result, in the limit of large wavenumbers considered here (i.e. 
for k — > oo), the third term in equation 1 1 3i can be neglected when 
compared with the first one. 

Finally, in the continuity equation < 141 . we note that all of the 
terms except the first one are multiplied by T or comparable terms 
[e.g. P' /{E + P)]. The latter is rather small and of the order ~ 
(c s /c) 2 , where c s and c are local sound speed and the speed of 
light, respectively. As a result, the leading terms we retain are the 



